Where is 99% of the condensation energy of Tl 2 Ba2CuO y coming from? 
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Anderson's interlayer tunneling model can account for up to 1% of the condensation energy of 
Tl2Ba2CuO y . Here we account for the remaining 99%. We predict an in-plane kinetic energy gain 
of 1 to 3 meV per planar oxygen in the optimally doped material when it goes superconducting. It 
is suggested that the effect may be most easily detected in underdoped dirty samples. 
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In conventional superconductors, the opening of the 
superconducting energy gap causes a reduction in optical 
absorption at low frequencies (below twice the gap value) 
U, and this missing spectral weight is transfered to the 
zero-frequency S— function that determines the penetra- 
tion depth H (so-called Ferrell-Glover-Tinkham (FGT) 
sum rule). It was proposed in Ref. 3 that in certain 
superconductors where pairing would arise from lower- 
ing of kinetic rather than of potential energy, the FGT 
sum rule would in appearance be violated, leading to 
the following novel phenomena: (1) a London penetra- 
tion depth shorter than would be expected from the low 
frequency change in optical absorption and/or the nor- 
mal state effective mass 0, and (2) a decrease in optical 
absorption at frequencies much higher than, and unre- 
lated to, the superconducting energy gap, possibly near 
infrared or visible. Recent experimental observations by 
Basov and coworkers in several high T c oxides § pro- 
vide striking evidence for the existence of this physics 
in high T c cuprates: for c-axis light polarization, the 
change in low frequency optical spectral weight accounts 
for just about 50% of the total spectral weight in the 
zero-frequency d— function, indicating that the remaining 
spectral weight comes from frequencies higher than the 
highest frequencies reached in the experiment, namely 
150 meV. Note that this energy is about 6 times the 
value of the gap, within which all the optical spectral 
weight change should be contained according to conven- 
tional BCS theory. || These findings are also consistent 
with earlier experiments by Fugol et al that report 
an anomalous decrease in optical absorption at frequen- 
cies in the visible range when high T c samples become 
superconducting. 

The well-known interlayer tunneling theory (ILT) 
of the cuprates predicts a lowering of the c-axis kinetic 
energy as the system goes superconducting, so that the 
considerations of Ref. 3 apply to that model for the case 
when light is polarized along the c direction. Realizing 
this, Chakravarty [0 applied the analysis of Ref. 3 to 
the particular case of c-axis conductivity and found (not 
surprisingly) that a lowering of kinetic energy in that 



direction would be consistent with an apparent violation 
of the FGT sum rule for c-axis light polarization, and 
hence the phenomena discussed above. For planar light 
polarization instead, no apparent sum rule violation is 
expected within the ILT model. 

However, within the ILT theory the kinetic energy low- 
ering should entirely account for the condensation energy 
of the superconductor In fact, in the analysis of c- 

axis transport of various high T c oxides where agreement 
with ILT theory is claimed |9|, e.g. La2- x Sr x CuOn, no 
distinction is made as to which frequency range in the op- 
tical absorption the 5-function spectral weight is coming 
from: the entire weight is attributed to kinetic energy 
lowering, which in turn implies that ILT theory would 
also be consistent with observation of no sum rule vi- 
olation However, the condensation energy of op- 
timally doped T^Ba^CuOy (T12201) is estimated from 
specific heat measurements to be 80 — lOO/zeV^ per CuOi 
unit 1 1 3 1 . Even in the extreme case where the entire 
weight in the superfluid condensate would come from ki- 
netic energy lowering the c-axis penetration depth re- 
quired would be less than 2^tm PJIc|l, which is signifi- 
cantly smaller than the measured penetration depth in 
that material, ~ 17/im | , , E5[ . In other words, the 
measured penetration depth implies that the kinetic en- 
ergy of c-axis motion in the superconducting state of 
Tl2Ba2CuO v accounts for less than 1% of its conden- 



sation energy |^fLJ,[15|. 

Still, we argue that Basov's and Fugol's experiments 
have clearly demonstrated that in high T c materials at 
least some of the superconducting condensation energy 
originates in kinetic energy. It is then natural to ask 
whether the rest might too. If so, the remaining conden- 
sation energy should arise from lowering of the in— plane 
kinetic energy. However, Basov and coworkers report || 
no apparent sum rule violation in the in-plane response 
within experimental accuracy so far, approximately 10% 
& 

This is however not surprising. We argue that the rea- 
son a sum rule violation has been detected in c-axis trans- 
port and not in in-plane transport so far is simply because 
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the materials are in the dirty limit for c-axis conduction 
and close to the clean limit for in-plane conduction jl7) . 
Note that the kinetic energy associated with the mea- 
sured in-plane penetration depth of T12201 of ~ 2000/1 is 
of order 20 — 60meV per C11O2 unit, of which the conden- 
sation energy is only a tiny fraction. Still, as we discuss 
below, we expect that an in-planc sum rule violation will 
be detectable. 

Within the model of hole superconductivity jl8j , pair- 
ing arises from lowering of kinetic energy in all directions 
when the system goes superconducting. Thus, the model 
predicts optical sum rule violation both in the in-plane 
and in the inter-plane response || . Because in-plane cou- 
plings are much larger than interplane ones, the bulk of 
the condensation energy in the model arises from kinetic 
energy lowering in the plane. Nevertheless, as we show 
in this paper, the model can readily account for the ob- 
served large sum rule violation in the c-direction and a 
much smaller sum rule violation in the plane, consistent 
with existing experimental results. It is predicted that 
with increased experimental accuracy an in-plane sum 
rule violation will be detected, most easily in the under- 
doped regime, which should readily account for the re- 
maining 99% of the condensation energy in Tl2Ba>2CuO y 
and other high T c cuprates. In fact, the model predicts 
a lowering of kinetic energy that is much larger than 
the condensation energy, that is partially offset by an 
increase in potential energy. 

Anderson remarks |l that for several high T c materials 
other than TliBa,2CuO v the condensation energy is com- 
pletely accounted for by the lowering of the c-axis kinetic 
energy as determined from the c-axis penetration depth, 
and concludes that for those cases "this agreement effec- 
tively rules out any intralayer theory of high T c " . This 
conclusion is logically flawed: the observations are en- 
tirely compatible with the existence of a much larger ki- 
netic energy lowering from in-plane motion together with 
an increase in potential energy upon pairing, as predicted 
by our model. 

The frequency-dependent conductivity in the super- 
conducting state is given by <j\ s (lS) = Dd(u>) + a^^uj), 
with a\s 3 (u>) the regular part. As discussed in Ref. 3, 
the superfluid weight in the zero-frequency d— function, 
D, is given by 



kinetic energy in the u direction as the system goes su- 
perconducting. The London penetration depth is given 
by A M = c/(8Z3 AI ) 1 / 2 , and an apparent sum rule viola- 
tion exists if it is shorter than what is obtained from just 
the low frequency missing area, SAf §f§, indicating the 
existence of 6A£ even if it is not detected directly. The 
degree of sum rule 'violation' can be characterized by the 
parameter 



V, 



" SA? + SA£ 



(3) 



and was found by Basov et al to be approximately 0.5 
in the c direction for several high T c materials including 
172201. 

For application to ILT theory the normal state kinetic 
energy in the c direction in Eq. (2) is assumed to be 
negligible, and Eq. (2) yields 
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with a and d in-plane and interplane lattice constants. 
The negative of the superconducting kinetic energy, < 
—T c > s , is equated to the condensation energy, e com j. 
For the case of Tl^Ba-iCuOy A c ~ 17/im, d = 11. 6A, 
a = 3.9A, eq. (4) yields t con d — 0.98//ey per CuOi unit, 
two orders of magnitude smaller than the value estimated 
from specific heat measurements Jj"3"| |. 

In the model of hole superconductivity, the kinetic en- 
ergy in the u direction is given by fig}] 
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with c ia a hole creation operator. The effective hopping 
in the u direction is given by 



ttf* = i„ + At M n 



(6) 



and it increases linearly with hole concentration n. To a 
very good approximation || the change in kinetic energy 
Eq. (2) is simply given by the expectation value of the 
correlated hopping term 



D„ = 5 A>; 



(i) 



with SAf the low frequency missing area, due to intra- 
band transitions, and 
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the missing area from much higher frequencies, related 
to interband transitions, v gives the volume of the unit 
cell, n indicates the direction of light polarization, and 
the right side of Eq. (2) gives the change in the carrier's 
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where the subindex a indicates that only the anomalous 
expectation values are to be included || . 

In the clean limit at zero temperature we have simply 



6A? 



2 2 
ire a 



2H 2 v 



E < -Ti > s 



(8) 



where the average of the single particle kinetic energy 
includes also the 'normal' expectation values of the At 
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term, yielding the renormalized single particle hopping 
Eq. (6). The degree of sum rule violation in that case is 
given by 



V, 



< rpAt > 



p < T* > s + < T** > s , a ■ 



(9) 



The full Hamiltonian for the model includes in addi- 
tion to the kinetic energy an on-site and nearest neighbor 
Coulomb repulsion: 

Hcoui = Hu + H v = n A n t\ + v X! nin i ( 10 ) 

i <i,j> 

and the condensation energy of the superconductor is 
given by 

tcond= [< T l >„ — < T* > s ]— < T At >^ Q — 

< H v > Sya - <H V > S:a = Ae t + e At + e v + e v (11) 

where the expectation values with subindex a indicate 
that only the anomalous contributions are to be included. 
In Eq. (11), a sum over the different directions is implicit 
for the terms that involve nearest neighbors. Quite gen- 
erally, all contributions to e conc i other than eAt are neg- 
ative: in the paired state, particles are closer together 
and experience larger Coulomb repulsion, and the single 
particle kinetic energy < T t > is optimal (most negative) 
in the normal state. A positive e conc i originates from a 
large negative contribution from < T At > S)a , he. a large 
kinetic energy lowering due to correlated hopping, which 
is to some extent compensated by the other terms in Eq. 
(11). In weak coupling, Eq. (11) is given by the usual 
form 



A 2 

tcond = -7rN(e F ) 



(12) 



with Ao the energy gap and N(cf) the density of states 
at the Fermi energy. 

In the model Eq. (5), each site denotes an oxygen atom 
in a Cu-0 plane. Consider the following parameters to 
describe the CuO planes of T12201 (case 1): U = 5eV, 
t = 0.19eV, At = 0.29eV. In the optimally doped case, 
with hole concentration n = 0.045 (per O atom), we ob- 
tain T c = 85K. Explicit calculation with the expressions 
given in Ref. 18 yields the results shown in Fig. [I]. The 
sum rule violation parameter Eq. (9) is V a = 3.4% for 
the optimally doped case, well below the current experi- 
mental accuracy. It increases up to V a = 6.5% in the un- 
derdoped regime. As we will show below however, the re- 
sulting sum rule violation in the c-direction can be much 
larger, consistent with observations. The condensation 
energy per O atom is 53/zeF in the optimally doped case, 
consistent with observations {e con d ~ per CuOi 

unit). Note that it results from a much larger kinetic en- 
ergy lowering due to correlated hopping, ca* = 1.2meV, 
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FIG. 1. Average in-plane kinetic energy per lattice site at 
zero temperature vs. hole concentration for case 1 parameters 
(given in the text). The solid, dot-dashed and dashed lines 
give the total kinetic energy, single particle and pair contri- 
butions respectively. The sum rule violation parameter V a in 
the clean limit (Eq. (9)) is also shown as a solid line (right 
scale). The critical temperature versus hole concentration, 
with maximum T = 85K, is indicated by a dotted line, (b) 
The various contributions to the condensation energy per O 
site (Eq. (11)) versus hole concentration. 



which is partially compensated by an increase in poten- 
tial energy and in kinetic energy from single particle hop- 
ping. 

The penetration depth in the ab plane in the clean 
limit is given by 
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< — T* > s + < -Tf* >.,„] 



(13) 



For the optimally doped case, we find < — T l a > s = 
-nmmeV and < -T Q At >= -0.610m e y, and with 
d = \1.QA Eq. (13) yields A a = 3697A However because 
there are other atoms between the CuO layers which are 
involved in the superconductivity it may be more ap- 
propiate to use d ~ 3 A in Eq. (13), which then yields 
A a = 1880A, consistent with experimental observations 
|pj| . Also, with the parameters chosen above, the density 
of states per CuC>2 unit is N(€f) = 1.25eV -1 , consistent 
with experimental estimates for T12201. 

Similarly, Fig. || shows results when a nearest neigh- 
bor repulsion is included (case 2). Here, U = 5eV, 
V = 0.65eV, t = 0.19eV, and At = 0.51eV is chosen to 
again yield T c = 85K for the optimally doped case. Here 
the kinetic energy lowering is larger, and yields a sum 
rule violation V a = 7.3% for the optimally doped case 
which increases to V a = 15% in the underdoped regime. 
We expect the actual parameters describing T12201 to be 
somewhere in the range spanned by these two examples. 
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FIG. 2. Same as Figure 1 for case 2 parameters, with 
non-zero nearest neighbor repulsion, given in the text. Note 
that the violation parameter V a as well as the kinetic energy 
lowering eAt are substantially larger than in Fig. 1. 



If the condensation energy was in fact somewhat larger 
than the current estimates the model would require a 
smaller value of the parameter t which would lead to an 
even larger effect. 

What about the c direction? The simplest assumption 
||l8f is that t and At are reduced by the same factor. If 
we take A c = 17 fim from experiment [flB) , and assume 
in Eq. (4) that the c-axis kinetic energy lowering due 
to At supplies 50% of the 6— function weight, as found 
by Basov (i.e. 50% sum rule violation), we obtain < 
T At >= 0.25/ieV (per O atom). For the parameters of 
case 1 above, we get an anisotropy factor 



< T At > 

< T At > 



2440 



(14) 



and the kinetic energy from single particle hopping scales 
approximately with the same factor This corre- 

sponds to an anisotropy in hopping amplitudes (or effec- 
tive masses) t a /t c = m*/m* ~ 50 The penetration 
depth in the clean limit in the c direction is given by 
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and yields \^ lean = 4.3/im. The observed penetration 
depth in the c direction is however a factor of 4 larger, 
which we attribute to the effect of disorder reducing the 
contribution of 5 At in Eq.(l) p§: 



SAi ~ 5Af ean x 7rA /(/i/T) 



(16) 



with an impurity scattering rate H/t c ~ 50 Ao, well in the 
dirty limit. In other words, the fact that the penetration 



depth we obtain from Eq. (15) is much smaller than 
the observed one, assuming this is due to the effect of 
disorder, is consistent with the observation that the sum 
rule violation in the c-direction is much larger than in the 
plane. Quite generally, under the assumption that low 
and high frequency missing areas have the same ratio in 
all directions in the clean limit Eqs. (3) and (16) yield 
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(17) 



hence a larger scattering rate in the c direction will nec- 
essarily be associated with a larger sum rule violation in 
that direction. 

As a further consistency check on the validity of this 
analysis we consider the normal states conductivities. 
For Ao = 13meV as obtained from our model, the above 
analysis implies Ti/t c = 5240cm - . For in-plane trans- 
port instead, Puchkov et al find J2^] a scattering rate 
1/t q = 560cm -1 , a factor of 10 smaller (which is con- 
sistent with a much smaller in-plane sum rule violation) , 
and a n (uj — » 0) ~ 5QOO(f2cm) -1 . Hence we expect in the 
c direction a conductivity that is smaller by a factor of ap- 
proximately T a /r c x m*/m* ~ 500, i.e. <r£ ~ 10(Ocm) -1 
and frequency-independent over several times Ao. This 
is consistent with observations. 

Note that the c-axis kinetic energy lowering for these 
parameters is only 0.5% of the total condensation energy, 
and only 0.02% of the kinetic energy lowering originating 
in in-plane motion. Results for the other set of parame- 
ters discussed above are similar. 

Will it be possible to observe a sum rule violation in 
the in-plane response? As shown above, the sum rule 
violation increases in the underdoped regime and hence 
it may be possible that the in-plane violation will be de- 
tected even if the system remains in the clean limit, es- 
pecially if the nearest neighbor repulsion is appreciable 
(case 2 above). Furthermore, by introducing disorder, 
for example by ion irradiation p2| , it may be possible to 
reduce the contribution of SAi to the superfluid density, 
moving the system towards the dirty limit, thus increas- 
ing the violation parameter. In this connection it should 
also be noted that experimental results in the c direction 
show a faster variation of the violation parameter with 
doping than found in Figs. 1 and 2 fl6|| . This however is 
simply explained by the fact that the system is moving 
towards the clean limit with increased doping, as evi- 
denced by the rapid increase in the optical conductivity, 
thus increasing the contribution of S Ai faster than given 
by the results in Figs. 1 and 2. Quantitative fits will be 
discussed elsewhere. 

In summary, the model of hole superconductivity pre- 
dicted an apparent optical sum rule violation || long 
before it was experimentally observed, which would be 
a manifestation of the reduction in effective mass and 
consequent lowering of kinetic energy that occurs upon 
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pairing. By assuming that T12201 is in the clean limit 
for in-plane conduction and in the dirty limit for the in- 
terplane response (which is consistent with observations 
), we found that the model can explain experimental 
observations. In the optimally doped case, the kinetic 
energy lowering from c-axis motion accounts for ~ 1% of 
the condensation energy and that of in-plane motion is 
20-50 times larger than the condensation energy. In par- 
ticular, the model can account for the 99% condensation 
energy that is missing in the alternative proposed expla- 
nation of the sum rule violation experiments, the ILT 
model J7|]. The fact that the predicted in-plane kinetic 
energy lowering is much larger than the condensation en- 
ergy is fortunate and should allow for its experimental 
detection. 
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